Spin dynamics and magnetic order in magnetically frustrated Tb2Sn207 



o 
o 

(N 



p. Dalmas de Reotier,^ A. Yaouanc,^ L. Keller,'^ A. Cervellino,^'Q B. Roessli,^ 
C. Baines,3 A. Forget,^ C. Vaju,i P.C.M. Gubbens,^ A. Amato,^ and P.J.C. King^ 

^ CEA/DSM/Departement de Recherche Fondamentale sur la Matiere Condensee, ,38054 Grenoble, France 
"^Laboratory for Neutron Scattering, ETH Zurich and Paul Scherrer Institute, 5232 Villigen-PSI, Switzerland 
'^Low temperature facilities group, Paul Scherrer Institute, 5232 Villigen-PSI, Switzerland 
^ CEA/DSM/Departement de Recherche sur I'Etat Condensee, 
les Atomes et les Molecules, 91191 Gif sur Yvette, France 
^Department of Radiation, Radionuclides & Reactors, 
Delft University of Technology, 2629 JB Delft, The Netherlands 
^Laboratory for Muon-Spin Spectroscopy, Paul Scherrer Institute, 5232 Villigen-PSI, Switzerland 
^ISIS Facility, Rutherford Appleton Laboratory, Ghilton, Didcot, 0X11 OQX, UK 

(Dated: February 6, 2008) 

We report a study of the geometrically frustrated magnetic material Tb2Sn2 07 by the positive 
muon spin relaxation technique. No signature of a static magnetically ordered state is detected while 
neutron magnetic reflections are observed in agreement with a published report. This is explained by 
the dynamical nature of the ground state of Tb2Sn207: the Tb''"'' magnetic moment characteristic 
fluctuation time is ~ 10~^° s. The strong effect of the magnetic field on the muon spin-lattice 
relaxation rate at low fields indicates a large field-induced increase of the magnetic density of states 
of the collective excitations at low energy. 

PACS numbers: 75.40.-s, 75.25.+Z, 76.75.+i 
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Magnetic materials with antiferromagnetically coupled 
spins located on triangular motifs exhibit geometrical 
magnetic frustration because their spatial arrangement is 
such that it prevents the simultaneous minimization of all 
the interaction energies 0|. The frustration, which leads 
to a highly degenerate ground state, forbids magnetic 
order to occur. Perturbations to the nearest-neighbor 
exchange interaction, such as exchange interactions ex- 
tending beyond nearest-neighbor magnetic atoms, dipole 
coupling or magnetic anisotropy, are believed to be re- 
sponsible for the magnetic order observed in some com- 
pounds 0. Typical examples are given by the spinel 
structure oxide LiMn204, the pyrochlore structure com- 
pounds Gd2Ti207, Er2Ti207, Tb2Sn207 and Gd2Sn207 
and the cuprate mineral Cu2Cl(OH)3. According to neu- 
tron diffraction (ND) measurements, long- and short- 
range orders coexist in LiMn204 j^, a partial order of 
the Gd'^+ magnetic moments is established at low tem- 
perature in Gd2Ti207 0] and a conventional magnetic 
order exists for the last three pyrochlore compounds 
[E IE Q' When looked for, persistent spin dynamics 
has always been found far below the magnetic ordering 
temperature. The possibility of such dynamics is con- 
ceivable for the partly ordered structure of Gd2Ti207 
and Cu2Cl(OH)3 but more exotic when all the mag- 
netic moments contribute to the magnetic structure as for 
Er2Ti207 0, Gd2Sn207 [lllk] and Tb2Sn207 U- A 
prerequisite for understanding the unanticipated behav- 
ior of these latter systems is a careful characterization of 
their dynamical properties. 

Here we show that positive muon spin relaxation (/iSR) 
and ND results in the ordered phase of Tb2Sn207 can be 
simultaneously accounted for only if the Tb'^^ moments 



are strongly dynamical. An independent and consistent 
time scale is obtained from a careful analysis of the neu- 
tron data. In addition, the initial strong and counter- 
intuitive increase of the muon relaxation rate when a 
magnetic field is applied indicates an increase of the den- 
sity of magnetic excitations at very low energy. 

Tb2Sn207 crystaUizes with the cubic space group 
Fd3>m. Rietveld refinements of powder x-ray and ND 
patterns yield the lattice constant a — 10.426 A and the 
free position parameter allowed by the space group for 
the 48/ site occupied by oxygen, x = 0.336 Mag- 
netic measurements point to a magnetic transition at 
0.87 K and to strong antiferromagnetic interactions as de- 
duced from the larg e and negative Curie- Weiss constant 
^cw = — 12K Powder ND indicates a structure 

with both ferromagnetic and antiferromagnetic compo- 
nents below Tgi = 1.3 (1) K where short-range magnetic 
correlations which are not liquid-like appear [6j . A steep 
increase of the Tb^+ magnetic moment /ixb and correla- 
tion length Lc is observed around T\r — 0.87 K where a 
peak is seen in the temperature dependence of the spe- 
cific heat Cp(T). ^xb = 5.9(1)^b with Lc ~ 19 nm at 
0.10 K. Therefore, even far below Tir, Lc is much shorter 
than usually observed in ordered magnetic structures. 

We present below (i) Cp{T) data recorded using a dy- 
namic adiabatic technique, (ii) ND measurements car- 
ried out at the cold neutron powder diffractometer DMC 
of the SINQ facility at the Paul Scherrer Institute and 
(iii) /^SR measurements done at the MuSR spectrome- 
ter of the ISIS facility (Rutherford Appleton Laboratory, 
Chilton, United Kingdom) and GPS and LTF spectrom- 
eters of the Swiss Muon Source (Paul Scherrer Institute, 
Villigen, Switzerland). 
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FIG. 1: (color online), (a). Low temperature dependence of 
the specific heat per mole of Tb measured for our Tb2Sn207 
sample. The sharp maximum is the signature of a magnetic 
transition occurring at — 0.88 K. The solid line is a fit to 
a model explained in the main text. (b). Magnetic powder 
diffraction pattern of Tb2Sn207 versus the scattering angle 
29 obtained from the subtraction of data recorded at 0.11 K 
and 1.23 K. Neutrons of wavelength 2.453 A were used. The 
solid lines show the best refinement and the difference spec- 
trum (bottom). From the pattern analysis we extract a Tb"^"*" 
magnetic moment fiTh = 5.4 (1) fiB which makes an angle of 
14.2° with the local (111) axis. A measurement (not shown) 
at 0.98 K gives /ixb = 1.6 (2) /ib and a Rietveld refinement of 
the diffraction pattern recorded at 100 K is consistent with 
space group Fd3m (lattice and oxygen position parameters 
equal to a = 10.427 A and a; = 0.337). 
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FIG. 2: (color online), (a). Two sets of /iSR spectra recorded 
respectively at 0.17 and 2.4 K. For clarity the latter ones have 
been vertically shifted by 0.05. The numbers written next to 
each spectrum correspond to the longitudinal magnetic field, 
expressed in tesla, for which the data have been recorded. 
The solid lines are fits to the stretched exponential relaxation 
function, P^'^{t) — exp[(— Azt)"], where a = 1 in zero-field, 
(b). Spin- lattice relaxation rate Az versus temperature, de- 
duced from spectra recorded using a cooling down sequence. 
The data were taken in zero-field, with the exception of some 
spectra above 4 K which were recorded in a longitudinal field 
Bcxt ~ 2 mT. All the corresponding spectra were fitted with 
a = 1. The solid line is deduced from the density of states 
introduced in the main text. 



We first present Cp{T) and powder ND measurements; 
see Fig.^ Relative to the published Cp{T) data 6], our 
sample displays a peak somewhat stronger in intensity 
at Tir = 0.88 K. The magnetic ND pattern at low tem- 
perature is in very reasonable agreement with published 
results The magnetic reflections are not resolution 
limited and their tails are Lorentzian-like rather than the 
usual Gaussian-like. Their shape has been fitted in detail 
assuming a size distribution of sharply defined spherical 
magnetic domains. A magnetic reflection is then rep- 
resented as the distribution-weighted sum of the peak 
profiles of domains of different sizes, each of them be- 
ing convoluted with the instrument resolution function 
Two types of distributions, namely the log-normal 
and gamma distributions, were used for the fit of the 
data shown in Fig. and gave similar results. For in- 
stance the average radii for the domains are 3.14 (10) 
and 2.95 (15) nm respectively for the two distributions. 
The volume-averaged domain diameter Dy = Lc can be 
deduced from the third and fourth moments M3 and Af4 
of the distribution using Dy = 3M4/(2M3) 0. Numeri- 
cally, wc find 19.8 and 19.1 nm for each of the considered 
distributions in excellent agreement with the value of Lc 
given above. 

Now we report on the /iSR data: see Refs. 0, 0| 



for an introduction to this technique. In Fig. |2K two 
zero-field spectra are presented, one recorded at a tem- 
perature T > Tsr and a second deep in the ordered state 
i.e. T ^ Tir. Unexpectedly, the two spectra are qualita- 
tively similar, i.e. no clear-cut signature of the magnetic 
transition below Ti^ is detected. We shall argue that this 
reflects the non-static character of the magnetic ground 
state of Tb2Sn207. 

We recall that the zero-field /iSR technique gives ac- 
cess to the longitudinal polarization function P^^^(t) and, 
in the magnetically ordered state of a powder sample, 
it is expected to be the weighted sum of the longitudi- 
nal and transverse components: P^^^{t) = [exp(— A^t) -|- 
2exp(— Ajcf) cos(7^(i3ioc)i)]/3. Xz and Xx are, respec- 
tively, the spin-lattice and spin-spin relaxation rates, 7^ 
is the muon gyromagnetic ratio (7^ = 851.615 Mrad s~^ 
T^^) and (Sioc) stands for the mean value of the local 
field at the muon site. Since all the measured spectra are 
exponential-like, this requires Ax — Xz and the oscilla- 
tory behavior to be absent. 

The first requirement is expected to be satisfied in 
the motional narrowing limit, i.e. when the dynamics 
of the Tb'^"'" moments is sufficiently fast; see for exam- 
ple Ref. The disappearance of the oscillations in 
the transverse component can have two origins: either 
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(Sloe) = or Te <C (7^-Bioc)~"^ where Tc is the character- 
istic fluctuation time of B\oc- We consider now these two 
possibiHties. 

The local field is built up from the dipole fields gen- 
erated by the Tb^+ magnetic moments. Taking into ac- 
count the magnetic structure, we have mapped the dipole 
field in the unit cell. The computed field is small enough 
to be consistent with the experimental result only in the 
neighborhood of site (0.212, 0.537, 0.463) and symmetry 
equivalent sites. This site is located at ~ 1.4 A from 
the closest oxygen atom neighbor. This corresponds to 
a much larger distance than the one usually adopted 
by the muon site in oxides (range 1.0-1.1 A) (see e.g. 
Refs. HO])- Therefore this possibility is considered 
very unlikely. 

Since the hypothesis (Bioc) = does not hold we 
now estimate a value of Tc from our data. We as- 
sume for simplicity a transverse stochastic field jump- 
ing between two opposite orientations j^^. Generaliza- 
tion of this model would not change qualitatively the 
result. Referring to the spontaneous fields measured for 
GdaTiaOy ^ and GdsSnaOr |i3, ^loc is estimated to be 
0.2 T. Since with our model Xz = T^^ioc'^c, we compute 
Tc = 8 X 10^^^ s from the measurements of \z at low 
temperature (Fig. ^p) . The whole /iSR analysis is con- 
sistent since the motional narrowing condition is fulfilled 

(7^SiocTc ~ 0.01 « 1). 

Now we have to understand the observation of mag- 
netic reflections in neutron scattering. The fact that 
these reflections can be indexed in the crystallographic 
structure of Tb2Sn207 implies that the scattering is elas- 
tic or nearly elastic. The energy resolution of the DMC 
diffractometer given by the energy spread of the incident 
neutrons is AE — 0.4 meV, a typical value for this kind 
of instrument. In other words, such a scattering exper- 
iment probes the magnetic structure with a time scale 
At = h/AE = 1.6 X 10-12 s. Since Tc > At the ND and 
/iSR results are compatible. Now, a time scale directly 
related to our diffraction measurements can be estimated. 
We use the relation E = H^k-^ /{2m„) defining the energy 
of a neutron of wavevector k and mass m„. From a mo- 
mentum width of 0.04 A~i deduced from the average in 
k space of the domain size distribution, a time scale of 
at most 2 x 10"^° s is obtained. It is rewarding that this 
value is of the same order of magnitude as Tc- 

Our result implies that the neutron scattering is not 
purely elastic, i.e. it occurs with a finite energy transfer. 
This dynamical order is consistent with the apparent re- 
duction of the Tb'^"'" moment deduced from the nuclear 
specific heat relative to the neutron determination Q. 

\z{T) is presented in Fig.[5j3. The rate is first temper- 
ature independent from 100 K down to '--^ 20 K and then 
starts to increase. This is consistent with the building 
up of pair-correlations at low temperatures [T^ . \z in- 
creases steadily as the sample is cooled. Remarkably, no 
sharp anomalies are detected at either Tgi- or Th-. \z{T) 
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FIG. 3: (color online). Field dependence of Xz measured 
at low fields for Tb2Sn2 07 using an experimental protocol 
explained in the main text. An amazing initial increase of Xz 
is observed at 3.2 K, i.e. far above Tsr, and down to 0.17 K. 
The shooting up of Xz at high field displayed for T = 2.4 K is 
also quite unexpected. The full line specifies the initial slope 
at 0.17 K which is 0.11 ^.s'^rsiT'^. 



presents an infiection point around Tif. Tsr seems to cor- 
respond to a change in the slope of Xz{T). Below 0.7 K 
Xz is observed to be only weakly temperature dependent 
down to the lowest measured temperature. That thermal 
behavior will be further discussed below. 

We have also recorded /iSR spectra under longitudinal 
fields (see Fig . Whereas the sequence of the mea- 
surements has no importance for T > Tgr, it matters for 
T < Tsr, suggesting some hysteresis. In the latter case, 
the present results were deduced from spectra obtained 
by first cooling the sample in zero field from T > Tgr and 
then using a field-increase sequence for the recording. In 
contrast to the zero-field spectra, a stretched exponential 
function (see caption of Fig. ^ with a < 1 is needed to 
fit these ones |23|- a (equal to 1 in zero-field) decreases 
smoothly with field to reach a value of respectively 0.8 
for T > and - 0.7 for T < T^r at 40 mT and then 
keeps this value up to 100 mT. Xz{Bcxt) is presented in 
Fig. El The observed initial increase at low fields is un- 
expected. The slope is large and of approximately the 
same value for the two temperatures probed below Ti,- 
and is still detectable in the paramagnetic phase up to 
3.2 K. Since Xz decreases at large fields for all the tem- 
peratures, an extremum takes place at intermediate fields 
for T < 3.2 K. Its position is located around 5 mT at 3.2 
K and goes up smoothly to reach a value ^ 40 mT at 
0.17 K, i.e. almost an increase of an order of magni- 
tude for more than a decade of temperature variation. 
The extremum being observed both in the ordered and 
paramagnetic states cannot originate from the response 
of magnetic domains or domain walls to -Boxt- 

The presence of such an extremum in Xz{Bcxt) bas 
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not been reported so far. The conventional Bloch- 
Wangsness-Redfield theory, which is expected to be vaUd 
at least for a temperature sufficiently large relative to Ti^, 
as at 3.2 K, predicts a monotonous decrease of Xz{BcKt) 
as JBext is increased . In addition, even if the low field 
part of the data is discarded, such a description breaks 
down because preliminary measurements, see Fig. |31 in- 
dicates Xz{Bext) to shoot up at large fields, that is, above 
1.0 T at 2.4 K. Note that this effect is not due to a change 
in a since a « 2/3 for all fields from 0.4 T upwards. 

We shall now provide a discussion for the Xz and Cp 
behaviors. As usual for geometrically frustrated mag- 
netic materials, Xz is only weakly temperature depen- 
dent in zero-field at low temperature. This is accounted 
for by a Raman scattering process involving two mag- 
netic excitations with a density of magnetic states char- 
acterized by an upturn at low energy and a small gap 
A proportional to the temperature i.e. A = ak^T. 
A simultaneous fit of both Xz{T) and Cp{T) in the or- 
dered state can be obtained using the density of states 
gmi^) — bf_ie^^^^ + bsh^^^"^ for e > A, where the first term 
gives rise to a temperature independent Xz and the sec- 
ond term has been chosen to reproduce our experimental 
T^/^ dependence in Cp{T). In Figs.^i and[2l3 are shown 
the results of this model with a = 0.02 and and bgh 
respectively equal to 0.057 meV^^^^ and 73 meV~^/^ per 
Tb atom. In addition we take the spectroscopic factor g 
= 2 and the anisotropy and exchange fields equal to Be = 
10 T and = 5 T, respectively (see Eq. 1 of Ref. 8). In 
fact, since, with the assumption a 1, Xz fe^/a^, any 
set of a and 6^ parameters satisfying bf^/a ~ 3 meV~^/^ 
will equally fit the data. The specific heat imposes a con- 
straint: if bfj, is too large the first term in the density of 
states gives a contribution to Cp{T). Therefore our com- 
bined set of data enforces 6^ ^ 0.06 meV~^/^ and hence 
a < 0.02. 

The increase of Xz at low fields is a definite signature 
of a field induced increase of the density of excitations at 
low energy. Because of the relatively large Zeeman energy 
on the system compared to A, the effect of the closure 
of the gap by the field is negligible. Indeed, at i?cxt — 10 
mT the Zeeman energy is 3.1 fieV which is large in com- 
parison to the computed value < 1 fj,eV for the gap at 
0.17 K. Experimentally, attributing the stronger relax- 
ation at low field to an increase of b^ and using the value 
of the slope from Fig. 13 for T < we estimate the 
relative change of &^ between and 10 mT to be 0.21. 

In conclusion, the magnetic order with a limited cor- 
relation length observed in Tb2Sn207 is dynamical in 
nature and characterized by a time ~ 10~^° s, much less 
than previously inferred |g] but in agreement with our 
analysis of the neutron results. This is about an order 
of magnitude longer than estimated for the analogous 
compound Tb2Ti207 which does not order [23,01 ■ O™' 
result implies that the magnetic scattering of the neu- 
trons is not purely elastic. Such a dynamical ground 



state is believed to be at work in the heavy fermion su- 
perconductor UPta jl^ . A sharp change in the dynamics 
at the temperature where a strong specific heat anomaly 
is detected has already been reported for the pyrochlore 
compound Yb2Ti207 which in contrast displays no mag- 
netic reflections 26]. Tb2Sn2 07 is a new case since both 
a specific heat peak and broadened magnetic reflections 
are detected at low temperature. Given the time range 
of the fluctuations, the moments are not as correlated 
as it was inferred previously 0. Neutron spin-echo ex- 
periments may provide an independent signature of the 
dynamical nature of the ground state of Tb2Sn207. This 
dynamics is related to the existence of an appreciable 
density of magnetic excitations at low energy. We have 
discovered that it can be increased by applying a small 
magnetic field. 
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